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"Number theory is not pure Mathematics. It is the Physics of the world of Numbers. " 

Alf van der Poorten. 

Abstract. In their 1995 paper, Jean-Benoit Bost and Alain Connes (BC) constructed a 
quantum dynamical system whose partition function is the Riemann zeta function C(/3), 
where /3 is an inverse temperature. We formulate Riemann hypothesis (RH) as a property 
of the low temperature Kubo-Martin-Schwinger (KMS) states of this theory. More precisely, 
■ the expectation value of the BC phase operator can be written as 

Ml) = *&1>l>-iW> 

where N q = Ylt=i Vh is the primorial number of order q and tpb a generalized Dedekind ip 
function depending on one real parameter b as 

"Ir *(«)-« n i^Ef- 

X2 pev, P \q 

Fix a large inverse temperature j3 > 2. The Riemann hypothesis is then shown to be equiv- 
alent to the inequality 

N g \MN q )\C(P ~ I) >e^loglogiV g , 
for q large enough. Under RH, extra formulas for high temperatures KMS states (1.5 < /3 < 
*0 ■ 2) are derived. 

1. Introduction 



X 



The Riemann Hypothesis (RH), that describes the non trivial zeroes of Riemann ( func- 
tion, is a Holy Grail of Mathematics [U |9]. Many formulations of RH may be found in the 



literature [5]. In this paper, we mainly refer to the Nicolas inequality [T3] 

N k 

> e 7 loglogiV fc , 



<p(N, 



where = Yli=iPi * s the primorial of order k, <p is the Euler totient function and 7 ~ 0.577 
is the Euler- Mascheroni constant. 
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The inequality ([I]) holds for all k > 1 if RH is true [lH Th. 2 (a)]. Conversely, if RH is false, 
the inequality holds for infinitely many k, and is violated for infinitely many k [TH Th. 2 (b)]. 

The Bost and Connes (BC) system is a quantum dynamical model whose partition function 
is the Riemann zeta function C(0) (with = the inverse temperature) and that exhibits 
a spontaneous symmetry breaking at the pole — 1. The behaviour of the system reminds 
that of a ferromagnet. At a temperature larger than the critical temperature T c w 10 3 -?T, 
the disorder of the spins dominates and the thermal equilibrium state is unique while for 
T < Tc the individual spins tend to align to each other, which in the three-dimensional 
space IR 3 of the ferromagnet yields a set of equilibrium thermal states parametrized by the 
symmetry group 50(3) of rotations in M. 3 . In Bost and Connes approach, the Hamiltonian 
operator is given by the logarithm of the number operator N as Hq = log N. This contrasts 
with the ordinary Hamiltonian H = N + | for the harmonic oscillator. It may be that the 
logarithmic Hamiltonian is realized in the field of quantum physiology, where the perception 
states are proportional to the logarithm of the excitation instead of being proportional to it. 
For a first account of this process, see [TU [12]. See also [17] for a model of a wave-packet 
dynamics based on a logarithmic energy spectrum. 

The purpose of the paper is two- fold. First, we rewrite the abstract mathematical language 
of BC theory with, as much as possible, concepts of quantum information processing (QIP). 
Second, we uncover a relationship between the statistical phase states (KMS states) of BC 
theory and Riemann hypothesis (RH). The Riemann hypothesis is shown to be equivalent 
to an inequality satisfied by the KMS states of the BC theory at a selected low temperature 
> 2. A special case (at — 3) of this inequality is ijj(N k )/(N k log log N k ) > e 7 /C(2) for 
k > 3, where ip(ri) = ip2 (n) is the Dedekind ip function [T3]. It plays a leading role for 
counting the number of maximal commuting sets of qudits [18]. Higher temperatures KMS 
states (1.5 < < 2) require RH to be analyzed with conventional tools such as the prime 
counting function ir(x). Very high temperatures KMS states (0 < < 1.5) are left as an 
open problem. 

2. On modular exponentiation: classical and quantum 

Let us start the parallel between the BC theory and QIP with the problem of finding the 
multiplicative order r in the modular equation 



(2) 



a x = 1 (mod q). 
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Given the integer a and a positive integer q with (a,q) = 1, the order r = ord,j(a) is 
the smallest positive integer x satisfying (J2J). It is straigthforward to establish the set of 
inequalities 

(3) OTd q {a)<X{q)< ( p{q)<q-l, 

where ip(q) is the Euler totient function [it counts the number of irreducible fractions | with 
p < q] and X(q) is the Carmichael lambda function [it is the smallest positive integer x 
satisfying §2§ for every integer a < q with (a,q) = 1] □. Both <p(q) and X(q) are universal 
exponents that satisfy (jSJ) for any integer a coprime to q. 

Order-finding is a task that can in principle be performed more efficiently on a quantum 
computer than on a classical computer. It is a building block of the Shor algorithm for 
factoring the integers [TU] . 

Let q be an m-bit integer and a be an integer smaller than q with (a, q) — 1. The quantum 
algorithm for order-finding is based on the following unitary operator acting on m-bits and 
such that for any n < q 

(4) \n) — > U a \n) = | an mod q) . 

Since a is coprime to q, it has an inverse modulo q and £7 a £^a = 1 (^o is unitary). Operator 
Z7 a is also multiplicative: £/ a [/& = UbU a , Va, b < q. 

The eigenvalues of the operator (j3J) are in the form exp(2i7r^) and the corresponding 
eigenvectors are given as a quantum Fourier transform 



1 ^— I —2mkj I • v 
(5) \ u k) = ^2^ ex v( : )Kmodg). 

The efficient implementation of the quantum Fourier transform on a quantum computer 
allows the efficient estimation of the period r of the function f(x) = a x mod q in (JSJ), 
and subsequently the efficient factorization of an integer (in polynomial time instead of the 
exponential time required on a classical computer). 



1 The recursive calculation of A(n) is as follows. For a positive integer k > 3, A(2 fc ) = 2 fe ~ 2 . For p 
an odd prime > 3 or < k < 2, X(p k ) = <p(p k ) = (p — l)p k ^ 1 . For q = Yl" =1 p ki a product of primes, 
X(q) = lcmfA^^Ji, where 1cm denotes the least common multiple of the numbers in the bracket. 
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3. On the KMS-states of Bost-Connes model 

Bost and Connes theory is a skilfull balance of concepts related to QIP and to quantum 
statistical physics [2]. The Bost-Connes C* algebra is the cross product [16] 

(6) Cq = C*(Q/Z) x N* 

of the C*-algebra C*(Q/Z) of additive/phase operators e$ acting on qudits |n) as 



(7) \n) -> eg \n) = exp(2ivr^) \n) , 6 = ^ e W z > 

by the semigroup N* of multiplicative operators \x a 



(8) |n) —> Ha \n) = \an mod q) , a G N*. 

The operator /i a is closely related to U a but acts on qudits |n). 

The BC system is the pair (Cq, <Jt) with entries the algebra Cq and a quantum statistical 
flow o t . 

In standard quantum statistical mechanics, given an observable Hermitian operator M 
and a Hamiltonian H , one has the evolution a t (M) versus time t as 

(9) a t (M) = exp(itH )Mexp(-itH ), 

and the expectation value of M at the inverse temperature (3 = is the unique Gibbs state 

(10) Gibbs(M) = Trace(Mexp(-/3iJo))/Trace(exp(-/3 J ff ). 

For the more general case of an algebra of observables A, the Gibbs state is replaced by the 
so-called Kubo-Martin-Schwinger (or KMS^ state). One introduces a 1-parameter group a t 
of automorphisms of A so that for any teK and x G A 

(11) a t (x) = exp(itH )x exp(—itH ). 

For the algebra A, the equilibrium state remains unique only if some conditions regarding 
the evolution of the operator x G A are satisfied: the KMS states tfip on A have to fulfill the 
KMS/3 condition with respect to a t , namely that for any x,y G A, there exists a bounded 
holomorphic function F x ^ y (z), < < f3 such that 

(12) F XtV (t) = M x °t(v)) and F XtV (t + i/3) = (/>p(<r t {y)x), Vt G R. 
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The KMS flow of BC system. The Hamiltonian on the BC system is defined from its 
action on the qudits \n) by 

(13) Hq \n) = Inn |n) . 

Using the relations exp(— ftH ) \n) = exp(— ftlnn) \n) = n _/3 \n), it follows that the partition 
function of the model at the inverse temperature ft is 

oo 

(14) Trace(exp(-/3tf )) = £>" " = C(/?), > 1- 

n=l 

Applying f TTTj) to the Hamiltonian ( 1T3"|) one obtains 

(15) cr t (n a ) = a %t n a and a t (e s ) = e s for any a G N*, 5 G Q/Z and t G R. 

It can be shown [21 [16] that in the high temperature regime < ft < 1 there is a unique 
KMS/? state of the quantum statistical system (Cq, c t ). It corresponds to the partition 
function (|14l) . In the low temperature regime ft > 1, the KMS/3 is no longer unique. The 
symmetry group G of the quantum statistical system is a subgroup of the automorphisms of 
at) which commutes with a t , i.e. 



(16) w o at = a t o w for any w E G and t G R, 

where o is the composition law. The symmetry group G permutes a family of extremal KMS^ 
states generating the possible states of the system after phase transition. The symmetry 
group of the algebra (u a , eg) endowed with the Hamiltonian H in (JTgj) is the Galois group 
G = Gal(Q c y cl /Q) of the cyclotomic extension Q c y cl of Q. The field Q c y cl is generated in 
Q by all the roots of unity. The Galois group G is aut(Q c y cl ). It is isomorphic to Z* 
At any temperature ft < oo, the restriction of the KMS^ state to C*(Q/Z) is given by 



(17) mi) =9-" n i 

and for ft = oo it is 



1 - p?- 1 



p^ 1 

pep Mi 



(18) 0oo(g) = M<?)Mg), 

where \x is the Mobius function. 

The KMS state f[T7j) may be easily rewritten as <fip(q) = N^Zi{q)^p-i{q), where N q 
n? =1 Pfc is a primorial number and ipb{<l) is defined in the abstract. 
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4. The Bost-Connes KMS states and Riemann hypothesis 

In sect ion f )4~Tj) one formulates Riemann hypothesis as a property of the KMS states (IT7|) 
with the inequality 



(19) 



e /3(?) 



N q \MN q )\ 



> when (3 > 2. 

j. In section (14.21) . new conditional 



logiogiv, CG0-1) 

The special case — 3 corresponds to our recent paper 
formulas are established. 

A concise account of numerical values at selected temperature are in table 1. 



q 


10 


10 2 


10 3 


10 4 


N q 


6.4 x 10 10 


4.7 x 10 219 


6.7 x 10 3392 


9.1 x 10 45336 


C2.1 


0.25 


0.093 


0.051 


0.031 


£3 


0.16 


0.018 


3.0 x 10~ 3 


6.1 x 10" 4 


eio 


0.25 


0.028 


4.9 x 10" 3 


1.0 x 10~ 3 


Table 1 . A check of KMS inequal 


ity (USD. 



4.1. Low temperature KMS states. To simplify the notation, we define the ratio Rb(n) := 
w log tog n ' m °ti va tion for this strange ratio will become apparent in the following appli- 
cation of Mertens formula [7] Th. 429]. 

Proposition 1. Pick a real number b > 1. For n going to oo, we have 

\im R b (N n ) 



C(bY 

Proof: By the definition of ipb{n) we can combine the Eulerian product for ((b) with the 
Mertens formula 

to obtain 

e 7 
C(6) 

Now the Prime Number Theorem [TJ Th. 6, Th. 420] shows that x ~ 9(x) for x large where 
9{x) stands for Chebyshev's first summatory function: 

p<x 
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This implies that, upon taking x = p n we have p n ~ @{p n ) = \og(N n ). The result follows. □ 
Proposition 2. Pick a real number b > 1. For n > 2 we have 



Ti 2 

n 2 > ip(n)il) h (n) > — 

Proof: The first inequality follows at once upon writing 

(p(n)ip b (n) 



n 2 



E[( 1 - VA 



p\n 

a product of finitely many terms < 1. Notice for the second inequality that 

(p{n)ip b {n) 



n 2 



- i /p b ) ^ ri( i - vp 6 



p\n p 

an infinite product that is the inverse of the Eulerian product for 

□ 

The next result is a far reaching consequence of Nicolas study of the small values of Euler 
ip function. 

Theorem 1. Under RH the ratio R b (N n ) is > ^ for n > 3. // RH is false, this is still 

true for infinitely many n. 

Proof: 

Follows by Proposition [21 combined with [HI Theorem 2]. □ 
In view of this result and of numerical experiments the natural conjecture is 

Conjecture 1. Fix a real number b > 1. For all n > 3 we have R b {N n ) > 
The main result of this note is the following. 

Theorem 2. Conjecture^ is equivalent to RH. 

Proof: If RH is true we refer to the first statement of Theorem [TJ If RH is false we consider 
the function 

Observing that log#(p n ) = logiV n , we see that g(p n ) < 1 is equivalent to Rb(N n ) > ^y. 
We need to check that for x large enough g(x) can be > 1 or equivalently logg(x) > 0. Using 
[21 Lemma 6.4], we obtain, upon writing 



- iogc(6) = - l /p b ) + E - Vp' 

p<x p>x 
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the bound \ogg(x) > \ogf(x) — b/x, where / is the function introduced in [TH Theorem 3], 
that is 



f(x) :=enog6{x)l[{l-l/p). 

We know by [Tl"| Theorem 3 (c)] that, if RH is false, there is a constant < c < 1 such 
that limsup x~ c f(x) > 0. 

Since for x large b/x << x~ c , the result follows. □ 

4.2. High temperature KMS states. To investigate high temperature KMS states at 
large dimensions n, we will need the following lemmas. Recall ir(x) the prime counting 
function 

7r(s)=5> 

p<x 

Lemma 1. Let < b < 1. We have 

p<x 

where J b (x) = z p|^. 

Proof: Write the above sum as a Stieltjes integral, and integrate by parts. □ 
Recall the Prime Number Theorem with rest from 1896 [7], 

ir(x) = Li(x) + 0(x exp(— c^/log x)), 
for some constant c > 0, as well as the Prime Number Theorem with rest under RH Q15J). 

tt(x) = Li(x) + 0(y/x\ogx). 

Plugging these estimates into the definition of J& we obtain 

Lemma 2. Let < b < 1. We have, unconditionally, 

" x Li(t)dt f x exp(-c v ^ogt)dt 



J b (x) 
and, under RH 



2 



t i + b v /, t b 



, Li(t)dt f x logtdt. 
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The preceding result shows that the assumption b > 1/2 is necessary to have asymptotic 
results to the precision 0(1), since we need the integral in the error term to converge. 

Lemma 3. Let 0.5 < b < 1. We have, under RH 

5Zi=^+w.(*)+o<D, 

p<x 

where 

Proof: Combine Lemma 1, Koch's paper [TBJ and Lemma 2. □ 
Theorem 3. Let 0.5 < b < 1. We have, under RH, 

Proof: Use lemma 3 and integrate lb by parts. □ 
If 0.5 < b < 1, denote by C& the sum of the convergent series 

a = !>>g(i - 1//) + i/A 

Denote by Ab(x) the partial sum 

A b (x) = J2^g(l-l/ P b ). 

p<x 

Denote by B b {x) the Bertrand integral B b (x) = f* t b\l gt - 

Theorem 4. Under RH and for 0.5 < b < 1, i/iere zs constant Kb > 0, stzc/i t/iat, for large 
n, we have 

ifo>(N n ) ~ K b logp n exjp(-B b (p n )). 

Proof: Note that the definition of Ab(x) entails 

A b (x) = C b + o{l)-J2W- 

p<x 

Combine Theorem 1 with the above and with Mertens formula 

JJ(1 - 1/p)- 1 ~ e 7 logx. 

p<x 
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The result follows upon putting x = p n . □ 
We now give a lower bound on ip b (N n ). 

Theorem 5. Under RH and for 0.5 < b < 1, for every e > 0, there is an N such that n > N 
yields 

iP b {n) > K b e\l - e)iV n log^ jV n exp(-£ 6 (p n )). 



Proof: First, note that 



<p(N n 



-MNn) =]J(l- 1/ P b 



N n 

P<Pn 

a divergent product bounded below by 

K b (l - e) log 2 N n exp(-B b (p n )), 

for n large enough, by Theorem HI 
Under RH Nicolas inequality 



>e^N n \og 2 N n 



holds for all n. The result follows. □ 

5. Conclusion and open problems 

In the present paper, we established the new criterion (fT5|) (theorem 2) for Riemann 
hypothesis from the low temperature KMS states of the Bost-Connes model. Our vision 
is related to but distinct from Connes vision of the zeroes of ( as an absorption spectrum 
[3]- In our approach the zeroes lurk in the backstage, and partial Eulerian products appear 
on the forefront. While the low temperatures KMS states ((3 > 2) lend themselves to an 
elegant analysis leading to infinitely many equivalence criterion with RH, high temperatures 
(2 > j3 > 3/2) forced us to assume RH to establish the asymptotic formulas in theorems 4 
and 5. 

Future work should focus on the physical understanding of the lower bound on KMS states 
(equivalent to RH) we established. Another direction of inquiry is the study of such a lower 
bound attached to a quantum statistical system whose partition function is the Dedekind 
zeta function [4j. As a finer point, we mention very high temperatures KMS states ( /3 < 3/2), 
which are challenging to analyze even under RH. 
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